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The  role  of  electronic  excitation  in  inelastic  collisions  between  ultracold  Ca  atoms  and  Ba+  ions, 
confined  in  a  hybrid  trap,  is  studied  for  the  first  time.  Unlike  previous  investigations,  this  system  is 
energetically  precluded  from  undergoing  inelastic  collisions  in  its  ground  state,  allowing  a  relatively 
simple  experimental  determination  and  interpretation  of  the  influence  of  electronic  excitation.  It  is  found 
that  while  the  electronic  state  of  the  ion  can  critically  influence  the  inelastic  collision  rate,  the  polar¬ 
izability  mismatch  of  the  neutral  atom  electronic  states  suppresses  short-range  collisions,  and  thus 
inelastic  processes,  involving  electronically  excited  neutral  atoms.  As  a  result  of  these  features,  it  is 
experimentally  demonstrated  that  it  is  possible  to  mitigate  inelastic  collision  loss  mechanisms  in  these 
systems,  marking  an  important  step  toward  long-lived  hybrid  atom-ion  devices. 
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Since  the  inception  of  laser-cooling,  a  chief  focus  of 
atomic  physics  has  been  the  development  of  techniques  for 
the  production  and  study  of  ultracold  matter — an  endeavor 
that  is  centered  on  facilitating  full  control  over  matter  at 
the  quantum  level.  This  work  has  been  quite  successful, 
enabling  many  long-sought-after  goals,  such  as  quantum 
degenerate  gases  [1,2],  quantum  simulation  [3],  quantum 
information  [4],  and  precision  measurement  of  fundamen¬ 
tal  physics  [5,6].  To  gain  this  control,  however,  most  ultra¬ 
cold  matter  production  techniques  rely  on  the  scattering  of 
a  large  number  of  photons  from  the  system  under  study, 
potentially  leading  to  a  large  degree  of  electronic  excita¬ 
tion.  Thus,  the  atom  or  molecule  being  cooled  can  be  in  a 
peculiar  state:  its  external  motion  may  be  characterized 
by  a  temperature  close  to  absolute  zero,  but  its  internal 
electronic  degree  of  freedom  may  be  characterized  by  a 
temperature  approaching  infinity.  While  this  non-thermal 
distribution  of  electronic  states  has  some  notable  important 
consequences  for  ultracold  atoms  [7],  most  of  its  effects, 
such  as  photochemical  reactions  [8-10],  are  largely 
ignored  as  they  occur  at  rates  that  have  relatively  little 
effect  on  the  system.  However,  as  the  field  now  moves 
towards  producing  more  complex  systems  at  ultracold 
temperatures,  e.g.,  molecules  [11]  and  hybrid  systems 
[12,13],  the  effect  of  these  light-assisted  processes  must 
be  reevaluated  as  their  rates  may  be  much  larger  due  to, 
among  other  things,  an  increased  density  of  accessible 
product  states  and  longer  range  interactions.  Thus,  there 
is  presently  a  need  to  better  understand  the  role  of  elec¬ 
tronic  excitation  in  inelastic  processes,  such  as  chemical 
reactions,  at  ultracold  temperatures  to  enable  the  next 
generation  of  experiments. 

Interestingly,  the  rapidly  emerging  field  of  hybrid 
atom-ion  systems  offers  a  unique  opportunity  to  study 
ultracold  collisions  [14-17].  Like  in  all-neutral  systems 


of  traditional  atomic  physics,  the  co-trapped  species  can 
collide  at  short-enough  range  for  reactions  to  proceed, 
but,  unlike  all-neutral  systems  it  is  possible  to  maintain  a 
product  of  the  reaction  in  the  trap,  since  ion  trap  depths  are 
large  relative  to  the  typical  exothermic  energy  release. 
Already,  such  hybrid  systems  have  been  used  to  measure 
several  important  ion-neutral  chemical  reactions  [13-16], 
culminating  in  the  recent  observation  of  molecular  ion 
reaction  products  [17]. 

Here,  we  use  a  hybrid  atom-ion  MOTION  trap  system 
[16,17]  to  study  chemical  reactions  between  ultracold 
40Ca  atoms  and  138Ba+  ions.  Unlike  all  previous  studies, 
inelastic  collisions  between  ground-state  atoms  and  ions 
are  energetically  forbidden  in  the  40Ca  +  138Ba+  system 
and  can  thus  only  occur  if  at  least  one  species  is  elec¬ 
tronically  excited.  As  a  result,  this  system  is  relatively 
simple  to  interpret  and  allows  a  straightforward  exami¬ 
nation  of  the  role  of  electronic  excitation.  In  what 
follows,  we  discuss  how  optical  coherence  and  trap  dy¬ 
namics  can  critically  affect  the  observed  atom-ion  chemi¬ 
cal  reaction  rate  when  electronic  excitation  is  involved, 
and  present  experimental  results  detailing  the  observed 
chemical  reaction  pathways  and  the  average  radiative 
association  branching  ratio  for  formation  of  CaBa+ 
molecular  ions.  We  also  present  relativistic  ab  initio 
molecular  structure  calculations,  which  further  illuminate 
the  excited  state  entrance  channels,  as  well  as  provide 
guidance  for  future  experiments.  We  additionally  present 
the  first  detailed  experimental  investigation  of  the  impor¬ 
tance  of  neutral  atom  excitation  and  conclude  with  a 
demonstration  that  it  is  possible,  with  reasonable  experi¬ 
mental  parameters,  to  avoid  the  large  chemical  reaction 
rates  seen  in  other  hybrid  atom-ion  systems,  marking  a 
crucial  step  towards  the  next  generation  of  hybrid  atom- 
ion  devices. 
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The  apparatus  used  in  this  work  is  similar  to  that 
described  in  Ref.  [16].  Briefly,  138Ba+  ions  (Nion  ~  1000 
ions)  are  loaded  into  a  linear  radio-frequency  (rf)  quadru- 
pole  trap  via  laser  ablation  of  a  BaCL  target  [  1 8]  and  cooled 
with  colinear  laser  beams.  In  the  same  region,  ultracold  Ca 
atoms  (A^mot  ~  3  X  106  atoms,  pMOX  ~  5  X  109  cm-3) 
are  produced  and  trapped  in  a  magneto-optical  trap 
(MOT).  (See  Fig.  1(a)  for  the  ion  and  atom  electronic 
energy  levels  and  cooling  laser  scheme,  denoted  by  arrows.) 
Two  nearly-orthogonal  cameras  allow  for  a  3-dimensional 
reconstruction  of  both  the  Ca  MOT  and  the  laser-cooled 
ion  cloud,  and  are  used  to  quantify  the  degree  of  overlap 
between  the  Ca  atoms  and  138Ba  ions,  as  shown  in 
Fig.  1(b).  Using  the  same  approach  as  Ref.  [16],  the  laser- 
induced  493  nm  fluorescence  monitors  the  138Ba+  ion 
population,  which  is  observed  to  decay  in  the  presence  of 
the  Ca  MOT  due  to  charge  exchange  chemical  reactions 
(Ca  +  Ba+  — > ►  Ca+  +  Ba,  CaBa+  +  y),  allowing  a  direct 
measurement  of  the  total  inelastic  collision  (charge  trans¬ 
fer)  rate  I  . 

In  this  experiment,  the  atom-ion  collision  energies  are 
primarily  set  by  the  coherent,  driven  motion  of  the  ion  due 
to  the  rf  trapping  potential  (i.e.,  the  micromotion)  since  the 
Ca  MOT  temperature  (T  —  4  mK)  and  the  Doppler-cooling 
temperature  limit  of  the  Ba+  ions  (T  —  0.5  mK)  are  both 

lower  than  the  typical,  position  dependent,  energy  of 

E 

micromotion  in  the  ion  cloud  (y21  ~  10  mK  —  10  K).  At 
these  low  collision  energies,  chemical  reaction  dynamics 
are  primarily  determined  by  the  relative  energies  of  the 
entrance  and  exit  reaction  channels,  shown  in  blue  (dotted) 
and  red  (dashed),  respectively,  in  Figs.  1(c)  and  1(d). 
Because  the  ionization  potential  of  Ca  (6.1  eV)  is 


significantly  larger  than  that  of  Ba  (5.2  eV),  charge  ex¬ 
change  reactions  between  the  Ca  4*5-  and  Ba+  6 2S-  or 
5 2D-  electronic  states  are  energetically  precluded.  As  a 
result,  inelastic  collision  between  Ba+  and  Ca  can  only 
occur  if  the  atom,  the  ion,  or  both  are  electronically  ex¬ 
cited.  Of  the  ten  energetically  allowed  entrance  channels 
populated  in  this  experiment,  the  six  involving  the  short¬ 
lived  Ca  P-states  do  not  contribute  to  the  observed  reaction 
rate.  Because  the  Ca  41 5'-,  4 lP-  and  5'P-states  have  very 
different  atomic  polarizabilties  [19]  (163.0  a.u.,  55.3  a.u., 
and  >1000  a.u.,  respectively),  during  the  collision  the 
strong  monopole  field  of  the  ion  alters  the  atomic  transition 
frequency  such  that  the  Ca  MOT  lasers  are  shifted  far  from 
resonance  for  atom-ion  separations  of  many  hundred  Bohr 
radii,  a0.  Thus,  though  a  collision  leading  to  reaction  may 
begin  with  a  Ca  atom  in  a  P-state,  by  the  time  the  atom  and 
ion  are  close  enough  to  react  (<  IOgoX  the  Ca  atom  has 
radiatively  decayed  from  the  P-state.  We  note  similar 
effects  are  discussed  in  Refs.  [14,20,21],  were  observed 
in  Ref.  [16],  and  are  confirmed  by  the  data  presented  here. 
This  effect  appears  general  and  it  is  likely  that  short-lived 
excitations  of  ultracold  atoms  can  be  ignored  in  future 
experiments  involving  reactions  of  ultracold  atoms  and 
ions.  (This  effect  is  particularly  beneficial  to  techniques 
for  sympathetic  cooling  of  molecular  ions  by  ultracold 
atoms  [22],  as  it  significantly  relaxes  the  experimental 
constraints.)  Interestingly,  a  similar  argument  cannot  be 
made  for  short-lived,  excited  states  of  ions.  Because  exci¬ 
tation  of  the  ion  only  changes  the  long-range  dispersion 
coefficient,  e.g.,  C6,  the  ion  cooling  lasers  are  not  shifted 
out  of  resonance  until  relatively  small  atom-ion  separa¬ 
tions.  Thus,  while  these  chemical  reactions  are  suppressed, 


Ba+  ion  cloud,  (c)  and  (d)  Molecular  potentials  of  the  fl  =  |,  |  excited  states  of  the  CaBa+  molecule,  respectively,  as  a  function  of 
interatomic  separation  R.  The  bottom  panels  show  higher  detail  for  potentials  near  20  X  103  cm-1  above  the  Ca+(25')  +  Ba(*D) 
dissociation  limit.  Entrance  and  exit  channels,  as  defined  in  the  text,  dissociate  to  blue  (dotted)  and  red  (dashed)  potentials, 
respectively. 
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it  is  nonetheless  possible  for  them  to  occur  between  atoms 
and  excited  state  ions. 

With  these  considerations,  it  is  clear  that  the  reactions 
proceed  from  the  energetically  allowed  entrance  channels 
involving  4*5-  or  3'D-state  Ca  atoms  and  6 2S-,  6 2P-, 
or  52D-state  Ba+  ions,  highlighted  in  green  in  Figs.  1(c) 
and  1(d).  While  the  charge  exchange  rate  constant  for  each 
of  the  four  distinct  pathways  depends  sensitively  on  the 
details  of  the  CaBa+  molecular  structure,  they  can  never¬ 
theless  be  probed  by  varying  the  excited  state  populations 
of  the  atoms  and  ions  while  recording  the  total  reaction 
rate.  Experimentally,  this  is  accomplished  by  varying  the 
intensities  of  the  atom  and  ion  cooling  lasers.  With  known 
experimental  parameters  it  is  then  possible  to  calculate 
the  relative  population  of  each  atom  and  ion  electronic 
state  (we  denote  the  population  fraction  of  the  specie  X  in 
the  electronic  state  i  by  its  diagonal  density  matrix  element 
pfi),  and  extract  the  contribution  of  each  reaction 
channel  to  the  total  measured  reaction  rate  given  by 
T  =  Pmot(°'v ) =  PMOT^Totai  =  P mot X V y'f >  where 
kjj  is  the  rate  constant  for  Ba+  in  electronic  state  i 
immersed  in  a  pure  sample  of  Ca  in  electronic  state  j. 

The  calculation  of  the  Ba+  electronic  state  populations 
is  complicated  by  optical  coherence  effects  and  the  quad- 
rupole  magnetic  field  of  the  MOTION  trap,  whose  direc¬ 
tion  and  strength  vary  in  a  non-trivial  way.  Therefore,  we 
determine  the  Ba+  electronic  state  populations  by  solving 
the  steady-state  quantum  Liouville  equation  for  an  eight- 
level  optical  Bloch  Hamiltonian  [23]  and  spatially  averag¬ 
ing  the  result  over  the  ion  cloud.  To  confirm  the  accuracy  of 
this  calculation,  we  directly  measure  the  Ba~62P|/2  state 
fraction  via  calibrated  fluorescence  and  ion  detection,  as 
shown  in  Fig.  2(a)  along  with  the  prediction  of  a  standard 
rate  equation  model  for  comparison.  Similarly,  a  first  prin¬ 
ciples  calculation  of  the  electronic  state  populations  of 
atoms  in  a  MOT  environment  is  a  well-known  open  prob¬ 
lem  and  various  phenomenological  models  have  been 
proposed  [24,25].  Therefore,  we  directly  measure  the  Ca 
4 2  P- state  excitation  fraction  to  determine  the  appropriate 


FIG.  2  (color  online).  Panel  (a):  experimentally  determined 
excitation  fraction  in  Ba+  versus  theoretical  predictions  of  the 
quantum  Liouville  (•)  and  rate  equation  (O)  treatment.  Panel 
(b):  experimentally  determined  excitation  fraction  in  Ca  versus 
predictions  of  a  4-level  rate  equation  model.  In  both  panels  the 
solid  line  has  a  unit  slope,  representing  ideal  agreement. 


model  for  calculation  of  the  Ca  population  fractions.  As 
shown  in  Fig.  2(b),  we  find  good  agreement  with  a  standard 
four-level  rate  equation  model  [10],  a  likely  consequence 
of  the  lack  of  hyperfine  structure  in  4<)Ca. 

With  these  tools  in  place,  we  vary  the  laser  intensities, 
calculate  the  resulting  excitation  fractions,  and  record  the 
total  inelastic  rate  constant.  The  resulting  measurements, 
plotted  as  a  function  of  ion  and  atom  excitation  fraction, 
are  shown  in  Figs.  3  and  4,  respectively. 

As  a  function  of  Ba+  excitation,  the  most  obvious 
feature  is  a  linear  trend  of  increasing  reactivity  with 
increasing  p®a  In  addition,  the  reaction  rate  constant 
diminishes  at  both  extremes  of  the  Ba+  S-state,  p®a,  and 
D-state,  pdd,  population  fraction.  Taken  together,  these 
features  indicate  that  the  primary  reaction  channel  involves 
the  Ba+  ion  P-state,  since  a  comparable  S-  or  D-state 
contribution  would  result  in  measurable  rate  constants  in 
the  limit  of  unit  occupation  of  the  S-  or  D-state.  Also 
shown  in  Fig.  3,  as  open  triangles,  are  several  measure¬ 
ments  at  low  excitation  fractions  (<  1%).  While  these 
points  agree  with  the  rest  of  the  data,  they  are  not  included 
in  the  final  analysis  to  avoid  possible  systematic  inaccur¬ 
acies  due  to  higher  secular  energies. 

As  shown  in  Fig.  4,  the  rate  constant  exhibits  no  detect¬ 
able  variation  over  the  experimentally  accessible  Ca  exci¬ 
tation  range.  Neither  the  Ca  52P-state  nor  the  41  P-statc  are 
expected  to  contribute  to  the  rate  constant  shown  here,  the 
former  due  to  a  negligible  population  and  the  latter  due 
the  aforementioned  radiative  suppression  effect.  As  a  result 
only  the  Ca  S-  and  D-states  contribute  to  the  observed 
reaction  rate;  thus,  we  conclude  that  the  lack  of  a  trend  in 
Fig.  4  is  an  indication  of  similar  reaction  rates  for  the 
atomic  S-  and  D-states  with  Ba+. 


FIG.  3  (color  online).  Measured  rate  constant  as  a  function 
of  Ba+  population  fractions.  Experimental  result  in  black.  Best 
fit  ({*„„ ksd, kpd, kdd}={ 4.2  ±  1.9,  <3,< 20,  <7} X  10- 10 cm3 s"1) 
to  the  four  entrance  channel  model  shown  as  open  red 
circles  (O).  For  these  measurements,  Ca  population  fractions 
were  fixed  at  pj/  ~  0.75,  p^p  ~  0.20,  pdd  ~  0.05. 
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FIG.  4  (color  online).  Measured  rate  constant  as  a  function  of 
Ca  population  fractions.  Experimental  result  in  black.  Best  fit 
({kps,  Kd,  kpd,  kdd}  =  {4.2  ±  1.9,  <  3,  <  20,  <  7}  X  lO"10  cm3  s'1) 
to  the  four  entrance  channel  model  shown  as  open  red 
circles  (O).  For  these  measurements,  the  Ba+  population  frac¬ 
tions  were  fixed  at  p®a  ~  0.4,  ~  0.1,  p dd  ~  0.5. 


To  extract  channel-specific  reaction  rate  constants  from 
these  data,  we  quantify  the  total  reaction  rate  constant 
in  terms  of  the  excited  state  fractions  of  both  Ba+  and 
Ca  as  £Total  =  p%p*^kps  +  p*?pcddksd  +  pddpddkdd  + 
P®pPddkpd.  By  accounting  for  the  conservation  of  proba¬ 
bility  (pss  =  1  —  ppp  —  pdd)  and  the  fact  that  Ca  atoms  in 
the  P-state  radiatively  decay  to  the  S-state  before  the  atom 
and  ion  are  close  enough  to  react  (p*fa  =  p[f  +  ppp),  the 
totality  of  the  data  in  Figs.  3  and  4  is  subject  to  a  multi¬ 
parameter  fit  along  the  manifold  {p®p,  pdd,  pdd}  to  deter¬ 
mine  the  contribution  of  each  reaction  pathway  to  the  total 
rate  constant.  From  this  fit,  the  reaction  rate  constant  of 
the  primary  entrance  channel,  Ba  (62/J|  /2)  +  Ca(4  .S’),  is 
found  to  be  kps  =  4.2(1. 9)  X  10-10  cm3s_1,  while  the 
best  fit  for  the  rate  constant  of  the  only  channel  out  of 
the  Ba+  ground  state,  Ba+(625'1^)  +  Ca(3'D),  is  con¬ 
strained  to  be  <  3  X  10“ 10  cm3  s_1.  Because  the  remain¬ 
ing  two  entrance  channels,  Ba+(62P^2)  +  Ca(31D)  and 
Ba+(52£>3/0)  +  Ca(21D),  require  simultaneous  excitation 
of  the  atom  and  ion,  our  experiment  is  relatively  insensi¬ 
tive  to  them  and  we  set  upper  limits  on  them  of  kpd  < 
2  X  10-9  cm3  s-1  and  kdd  <7X  10“ 10  cm3  s_1.  To  show 
the  quality  of  agreement  between  the  best  fit  and  the 
results,  we  plot  a  data  point  generated  from  the  fit  corre¬ 
sponding  to  each  experimental  condition  (open  circles)  in 
Figs.  3  and  4.  In  addition  to  this  fit,  using  the  method 
described  in  Ref.  [16],  we  measure  the  molecular  ion 
product  branching  ratio  to  be  ~0. 1 .  Future  measurements, 
using  time-of-flight  detection  as  in  Ref.  [26],  are  planned  to 
further  constrain  the  branching  ratio. 

To  illuminate  the  charge  exchange  mechanisms  respon¬ 
sible  for  the  individual  reaction  pathways,  we  performed 


ab  initio  calculations  of  the  CaBa+  potentials,  with  total 
body-fixed  angular  momentum  of  fl  =  (^,|),  using 
a  relativistic  multi-reference  restricted  active  space 
configuration-interaction  method  [27],  Spin-orbit  effects 
are  large  in  CaBa+  and  a  relativistic  calculation  is  required. 
The  resulting  Born-Oppenheimer  potentials  for  the  ground 
and  first-excited  state  of  CaBa+  are  shown  in  Figs.  1(c) 
and  1(d),  and  assigned  by  their  atomic  dissociation  limit. 
The  bottom  panels  of  Figs.  1(c)  and  1(d)  show  a  higher 
resolution  view  of  the  closely  spaced  potentials  above 
20  X  103  cm-1.  Both  panels  reveal  that  there  are  strong 
interactions  and  avoided  crossings  between  neighboring 
excited  potentials.  Potentials  dissociating  to  energies 
>27  X  103  cm-1  are  not  shown,  as  the  exceedingly 
high  level  density  makes  it  impossible  to  resolve  the 
avoided  crossings  with  our  calculation. 

In  principle,  a  coupled-channels  calculation  that 
includes  the  mixings  between  these  potentials  due  to  radial 
non-adiabatic  coupling  and  spontaneous  emission,  like  that 
performed  in  Ref.  [16],  can  be  used  to  calculate  the  non- 
radiative,  radiative,  and  radiative  association  charge  ex¬ 
change  reaction  rate  constants  for  each  entrance  channel; 
however,  due  to  the  high  level  density  in  the  excited 
CaBa+  molecule,  this  calculation  is  extremely  technically 
demanding  and  falls  outside  the  scope  of  this  work. 
Nonetheless,  the  basic  mechanisms  of  each  charge  ex¬ 
change  pathway  can  be  qualitatively  inferred  from  the 
calculated  structure  of  its  entrance  channel  molecular  po¬ 
tential.  For  example,  the  experimentally  dominant  entrance 
channel,  Ca(4 1 S)  +  Ba+ (62/J]/2),  exhibits  strong  interac¬ 
tions  with  several  potentials  as  well  as  a  narrowly-avoided 
crossing  with  a  molecular  potential  dissociating  to  the 
Ca+(425l)  +  Ba(63P)  exit  channel.  The  second  most 
experimentally  reactive  entrance  channel,  Ca(3 1 D)  + 
Ba -  (62.Sj/2),  also  shows  several  strong  avoided  crossings 
with  the  Ca+(32D)  +  Ba(615l)  exit  channel.  Strong, 
narrowly-avoided  crossings  facilitate  non-adiabatic,  non- 
radiative  charge  transfer,  the  likely  mechanism  for  the 
aforementioned  channels.  By  similar  reasoning,  it  is 
expected  that  the  Ca(3 1 D)  +  Ba  (52D)  entrance  channel 
would  be  relatively  inefficient  for  charge  exchange  reac¬ 
tions  since  low  energy  scattering  events  do  not  experience 
an  avoided  crossing  nor  do  they  couple  well  to  states 
accessible  via  spontaneous  emission  due  to  poor  Franck- 
Condon  overlap.  The  remaining  Ca(3 1 D)  +  Ba+(6 2Pl/2) 
entrance  channel  resides  in  a  host  of  densely  packed  levels, 
which  prevents  our  ab  initio  calculation  from  providing 
reliable  results  and  suggests  that  non-radiative  and  radia¬ 
tive  chemical  reactions  could  proceed  very  quickly  from 
these  states,  a  result  that  is  consistent  with  the  limit  set  for 
this  channel  by  our  measurement. 

Finally,  despite  the  fact  that  the  rate  constants  measured 
here  are  large  by  atomic  physics  standards,  roughly  a  third 
of  the  rate  attained  when  assuming  all  collisions  that  cross 
the  angular  momentum  barrier  have  a  unit  probability  of 
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being  inelastic  (i.e.,  the  Langevin  rate),  with  proper  engi¬ 
neering  of  the  atomic  and  ionic  excited  state  populations, 
their  effects  can  be  mitigated.  By  using  sufficient  repump¬ 
ing  laser  powers  to  maintain  a  negligible  Ca  D-state  frac¬ 
tion  and  using  the  minimum  ion  laser  cooling  power 
necessary  to  achieve  Doppler-limited  cooling,  we  have 
observed  Ba+  ion  lifetimes  in  excess  of  several  minutes 
in  the  presence  of  the  MOT. 

In  summary,  we  have  experimentally  studied  chemical 
reactions  in  a  hybrid  atom-ion  system,  where  only  elec¬ 
tronically  excited  states  are  permitted  to  react.  The 
observed  reaction  dynamics  are  qualitatively  different 
than  what  is  observed  in  systems  where  ground-state  reac¬ 
tions  are  equally  important  [16].  The  measured  reaction 
rate  constant  is  found  to  be  an  appreciable  fraction  of  the 
Langevin  rate,  similar  to  results  found  in  other  experiments 
involving  excitation  [14,17,28],  and  is  most-likely  due 
to  non-radiative  charge  transfer  occurring  at  narrowly- 
avoided  crossings  in  at  least  two  of  the  entrance  channels. 
We  have  also  demonstrated  that,  with  a  judicious  choice  of 
atomic  species  and  degree  of  electronic  excitation,  it  is 
possible  to  engineer  experiments  where  the  cold  reactive 
species  can  co-exist  for  several  minutes.  From  these  results 
it  is  clear  that  a  new  generation  of  hybrid-atom  ion  devices, 
capable  of  the  long  coherence  times  needed  for  proposed 
hybrid  atom-ion  systems  [22,29-33]  is  in  reach. 
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